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thermal expansion of the spin-ladder system (C5Hi 2 N) 2 CuBr4 
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We present high-resolution measurements of the c*-axis thermal expansion and magnetostriction 
of piperidinium copper bromide (C5Hi2N)2CuBr4. The experimental data at low temperatures is 
well accounted for by a two-leg spin-ladder Hamiltonian. The thermal expansion shows a complex 
behavior with various sign changes and approaches a 1/s/T divergence at the critical fields. All low- 
temperature features are semi-quantitatively explained within a free fermion model; full quantitative 
agreement is obtained with Quantum Monte Carlo simulations. 



In recent years there has been an increasing inter- 
est in spin compounds exhibiting magnetic field-induced 
quantum phase transitions. Examples are a Bosc- 
Einstein condensation of magnons |l|, |2| studied in 
three-dimensional (3D) systems like coupled spin dimcrs 
[3j, |j, |5| or arrays of coupled spin-1 chains [y, |7| or spin- 
ladders [a |9j. These systems share a zero-field ground 
state with a spin gap. In a magnetic field two quantum 
phase transitions are present; at H c \ the gap closes and 
at H C 2 the fully field-polarized state is reached. An in- 
teresting scenario arises in the limit of very weakly cou- 
pled ladders or chains where an extended temperature 
regime controlled by ID physics is expected. The low- 
dimensionality makes the investigation of quantum crit- 
ical properties particularly exciting, and the thermal ex- 
pansion a is especially suited for this purpose. In fact, it 
has been shown that the Griineisen parameter T — a/C, 
where C is the specific heat, necessarily diverges at a 
quantum phase transition [101 ]. In addition, a shows a 
sign change whose location in the phase diagram indi- 
cates accumulation of entropy [111 ]. For D < 1/u, where 
v is the correlation length exponent, even a will diverge 
at criticality. For a spin-ladder, v = 1/2, the thermal 
expansion at the critical fields is expected to behave as 
a ~ 1/vT |ipj . Experimentally, this is largely unex- 
plored due to the lack of suitable model materials. 

In this letter we present a study of the thermal expan- 
sion and the magnetostriction of single crystalline piperi- 
dinium copper bromide (CsHj^N^CuBi^ whose mag- 
netic subsystem is a very good realization of a two-leg 
spin ladder [8j] with Hamiltonian 

N 

"H = 2_^ [-/xS^iS^a + J\\ (Si.iSi+i,i + Si,2Si+i,2) 

z=l 

-cjUbH (SI, + S? 2 )] . (1) 

The critical fields H cl = 6.8 T and H c2 = 13.9 T im- 
ply for the exchange couplings the values J±/ks = 12.9 



K and Jn/kB = 3.6 K 12]. The inter-ladder couplings 



are very weak since no indications of three-dimensional 
magnetic ordering are present down to T ~ 100 mK [151 ] . 



On approaching each of the critical fields, we find highly 
anomalous temperature dependencies of a. For 0.3 K < 
T < 2 K, a(T) approaches 1/VT divergences with op- 
posite signs for H c i and H c i in agreement with the ex- 
pected quantum critical behavior. Away from the critical 
fields, a shows a rather complex structure with various 
sign changes, essentially antisymmetric with respect to 
(H c i + H C 2)/2. All these low-temperature features are 
reproduced semi-quantitatively within a model of free 
fcrmions, although the latter is known to be valid only 
in the vicinity of J? c i/c2- 

(C5Hi2N)2CuBr4 crystallizes in a monoclinic struc- 
ture |lj] . The legs of the spin ladders are oriented along 
the a axis and the rungs roughly (w 20°) along c*. The 
single crystals used in this study have been grown from 
solution |l5(. We present high-resolution measurements 
of the thermal expansion a{T) = 1/L ■ dAL{T)/dT, 
magnetostriction s(H) = AL(H)/L and its field deriva- 
tive A = de(H)/dH. Here, Lq is the sample length along 
c*; AL(T) and AL(H) denote the temperature- and field- 
induced change at constant H and T, respectively. The 
measurements have been performed on a home-built ca- 
pacitance dilatometer in magnetic fields H\\c* up to 17 T 
for 0.3 K < T < 10 K. 

Fig. Q] displays e(H) measured at various constant tem- 
peratures. At T — 335 mK, e(H) is field independent up 
to about 6 T, then it continuously increases until it sat- 
urates for H > 15 T. This behavior strongly resembles 
that of the low-temperature magnetization |8( . The in- 
set of Fig. Q] shows the magnetostriction coefficient \{H) 
with two pronounced peaks at 335 mK, whose magni- 
tude shrink with increasing temperature until around 2 K 
both peaks merge into one broad plateau which further 
broadens towards higher T . In order to identify the crit- 
ical fields specified above we used linear interpolations of 
the peak positions of X(H) towards T = K. 

In Fig.[5{a) we summarize the thermal expansion data 
for H < H c \. In zero field, a(T) has a pronounced peak 
around 5 K followed by a strong increase above about 
15 K. The latter is due to the usual thermal expansion 
of phononic origin, while the 5 K peak is well described 




Magnetic field (T) 

FIG. 1: (a) Magnetostriction e = AL(H)/L and (b) A = 
de/dH measured along c* at various temperatures. The solid 
line in (a) is a fit of e at 335 mK within a free fermion model. 



by a Schottky formula indicating the thermal occupation 
of triplet excitations, see discussion below. These de- 
generate triplets are split upon increasing magnetic field 
which is reflected in a broadening and slight shift of the 
Schottky peak towards lower T. For H > 5.5 T, the mag- 
nitude of the peak starts to grow until around H c i a con- 
tinuous increase of a(T) down to our lowest T ~ 0.3 K 
is observed. For H > H c i, this low-T increase rapidly 
weakens and a becomes even negative for H — 7.2 T, 
see Fig. [2jb). With further increasing field this sign 
change shifts towards higher T and a clear minimum- 
maximum structure of a(T) becomes visible, whose am- 
plitude continuously decreases until a comparatively fea- 
tureless a(T) curve is reached at the aforementioned sym- 
metry field (H cl +H c2 )/2 = 10.4 T. On approaching H c2 
we observe the same systematics of a(T) but with in- 
verted signs, Fig.[5fc); in particular, a negative Schottky 
peak develops above H c2 , see Fig. (2fd). 

We now turn to the theoretical discussion of the mag- 
netic contributions to the thermal expansion and magne- 
tostriction arising from the Hamiltonian (p}. As J\\/J± ~ 
1/4 for (C5Hi 2 N) 2 CuBr4 we can resort to the strong cou- 
pling limit Jii <C Jj_ where the physics emerges in a trans- 
parent way. In this limit, the Hamiltonian |T]) represents 
a sum of weakly interacting dimers in ID. For k B T ^$> Ju , 
the magnetic contribution to a is then attributed to the 
thermal occupation of single-dimer states described by 
the Schottky formula 



CtSch 



ie P(h+J ± ) £ g _ e 2/3h) + J x (J + e Ph + e 2/3fc)] 

k B T 2 (l + ef ih + e 2 ^ + e^ h+J ^) 2 

(2) 
where 1//3 = k B T and 7 = y-^g^- with the volume per 

dimer V D ~ 859 A 3 15]. A fit for H = (dashed line in 
Fig-EI a )) yields 7 ~ 11.7 x 10~ 5 [16] implying a uniaxial 
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FIG. 2: Thermal expansion a = 1/L ■ dL/dT measured along 
the c* direction in various magnetic fields (a) below H c \, (b,c) 
between H cl and H C 2, and (d) above H C 2- Note the logarith- 
mic temperature scale. The dashed line in (a) is a fit with the 
Schottky formula Eq. ([2]). The arrows signal the order of the 
a(T) curves with increasing field. 



pressure, p c *||c*, dependence d\nJ±/dp c * ~ 55%/GPa. 
In lowest order, the intra-chain coupling just leads to 
a mean-field shift of the effective magnetic field, h = 
g[i B H — J\\M/{c)Hb), where M is the magnetization of a 
single dimer. The Schottky formula @ accounts well for 
the field dependence of the 5 K peak (not shown) . In par- 
ticular, it identifies for any given temperature a unique 
magnetic field where the thermal expansion vanishes and 
changes sign, see Fig. [3] Upon decreasing temperature, 
the positions of vanishing a in the (H, T) phase diagram 
shift to lower fields towards the location of the singlet- 
triplet transition, h ~ Jj_. 

However, when temperature is of order k B T ~ J|| the 
behavior changes qualitatively and the Schottky formula 
ceases to be valid. In contrast to a simple singlet-triplet 
level crossing of isolated dimers, the kinetic energy of the 
triplet excitations gives rise to an extended gapless phase 
at zero temperature terminated by two quantum critical 



points at H, 



cl/c2- 



In the low-T limit k B T < J±, the 



Hamiltonian (Q} can be mapped onto an effective XXZ 



spin-1/2 chain [19|, l20|, |2l| that is described by inter- 



acting tight-binding Jordan- Wigner fermions with band- 
width Jii and chemical potential tuned by magnetic field, 
M = 9^bH — J±+0(Ji\). Here, the magnetic field and the 
pressure dependence J±{p) only enter via the chemical 
potential \i which implies essentially the same behavior 
for the derivatives of the free energy with respect to H 
and p. This not only explains the close correspondence 
between magnetization [8| and e(H) of Fig. Q] in the low- 
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FIG. 3: Regions of positive and negative thermal expansion 
a along the c* axis of (CsHi2N)2CuBr4. The triangles and 
circles denote the positions of vanishing a determined from 
an interpolation along the T and H axis, respectively. In 
the shaded area a is negative. The dashed line indicates the 
zeros of a as predicted by the Schottky formula ([J}. The 
solid lines show zeros of a derived from the free-fermion fitting 
formula Q , which for T — start at the two quantum critical 
points located at H c \ = 6.8 T and H C 2 — 13.9 T and at 
(H c i + H C 2)/2. Note that a crossing point of lines of zeros as 
suggested by Eq. Q is not realized in the experiment as the 
regions of negative thermal expansion are connected. 



T limit but it also provides an intuitive interpretation of 
the various sign changes of the thermal expansion a. 

The sign change of thermal expansion, a ex dS/ dp, 
reveals the locations of extrema of entropy, S, in the 
phase diagram that originate from the proliferation of 
low-energy fluctuations close to a quantum critical point 
ll| . In the present context, the close vicinity of two 
quantum critical points H c \i c i gives rise to a rich struc- 
ture in a, see Fig. [21 whose sign changes we summa- 
rized in Fig. [3] At elevated temperatures, thermody- 
namics is not able to resolve the distance between the 
critical points, viz. the bandwidth Jn of triplet excita- 
tions. For temperatures Jn < k B T <C J±, this is reflected 
in a broad single peak in the entropy as a function of 
magnetic field S(H), or, equivalently, as a function of 
the chemical potential /x, indicating the singlet-triplet 
level crossing. This peak implies a single sign change 
of a ex dS/dJ± ex dS/d/i in agreement with the predic- 
tion of the Schottky formula ([2j). At k B T ~ Jn, however, 
thermodynamics starts to resolve the triplet bandwidth 
which is reflected in a splitting of the entropy peak into 
two separate maxima whose positions approach the crit- 



ical points H cl / c2 for T — > K. The bifurcation of the 
entropy peak is the origin of the rich structure of a. The 
two maxima and the enclosed minimum of S(H) result 
in three consecutive sign changes of a{H) for k B T < Jn 
as shown in Fig. [3J 

Close to criticality, H rs H c i/ c2 and k B T <C J\\, 
the interactions among Jordan- Wigner fermions can be 
neglected and the critical model reduces to free non- 
relativistic fermions. The resulting singular part of the 
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FIG. 4: Panels (a-c) compare representative experimental 
data (symbols) of a(T) to theoretical curves (lines) obtained 
within the free fermion model (|4]). All these curves are fully 
determined by H cl / c2 and a factor 7 obtained from a fit to 
e{H) data at T = 335 mK (see Fig. [TJ. Panel (d) shows the 
thermal expansion near H cl / c2 on a double-logarithmic scale 
together with the limiting 1/yT behavior (dashed line). 



free energy has the scaling form 






\gH B (H - H c 



|3/2, 



k R T 



g^B (H - H c 



(3) 



where n = 1,2 label the two distinct critical points and 
the scaling function 



was specified in [17 



that the critical fields H, 



l/c2 



Assuming 
are smooth functions of 



pressure, it follows from ((3]) that s(H) oc dF cr , n /dp <x 
dF crt7l /dH cn approaches the same characteristic square 
root behavior as the magnetization in the zero tempera- 
ture limit, e(H) ~ y/\H — H cn \; this is already apparent 
in the low-T data in Fig. [TJ For the critical thermal ex- 
pansion the previously announced divergence a ~ 1/yT 
is found for \H — H cn \ ^C T. Note that a divergent ther- 
mal expansion at a quantum critical point is not pro- 
hibited by any fundamental principle. However, $S§ also 
predicts a divergent critical contribution to the compress- 
ibility that might result in a preemptive first-order tran- 
sition of the elastic system [181 ] . 

In order to capture the complex behavior of a(T,H) 
for fc^T < J u we employ a simple fitting formula 



OLFF 



"7 



dk 



O-e/c) 



% 2ir Ak B T 2 cosh 2 ((/i - e k )/(2k B T)) 



(4) 



where p, — gpsi^H — H c \ — H C 2)/2, the dispersion 
e k — 2~gHB{Hc2 — H c i)cos(k), and 7 measures the pres- 
sure dependence of the singlet-triplet crossing field, 7 = 
g/j, B (d(H c i +H c2 )/dp)/(2V D ). Formula gj) follows from 
a free fermion model with a pressure dependent chemical 
potential fi and dispersion e^; note that it reproduces the 
qualitatively correct 1/yT behavior close to the critical 
fields. The fit of e{H) at 335 mK, see Fig. [JJ identi- 
fies 7 ~ 12 x 10~ 5 with the saturation value of e(-ff) at 
large fields, in good agreement with the value obtained 
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FIG. 5: Anisotropic magnetostriction coefficients A (symbols) 
measured along the c* , a, and b axis. The lines are Quantum 
Monte Carlo fits using the experimental values of J± and 
J|l and two additional parameters for each direction, which 
measure the uniaxial pressure dependencies of J± and Jii [22| . 



from the Schottky fit of a(T). Having 7 fixed, Eq. Q 
strikingly reproduces the complex behavior of a(T) in the 
low-T limit for all magnetic fields without further adjust- 
ment of any parameter, see Fig. [4] In particular, Eq. (|4]) 
confirms the three consecutive sign changes of a at low- 
est T as discussed above. The fit also demonstrates in 
Fig. |U[a) that in the measured temperature range the off- 
critical curve at H = 6.6 T is even larger than the critical 
one at H c \ = 6.8 T, but is expected to cross the latter at 
lower T. In Fig. EJd) we show the a(T) curves that are 
closest to the two critical fields on a double-logarithmic 
scale. Whereas a at H = 6.8 T is critical, the one at 
H = 14 T is slightly off-critical. The dashed line indi- 
cates the limiting 1/VT behavior. The deviation of the 
latter from the critical H = 6.8 T fitting curve (solid line) 
is attributed to certain corrections to scaling that are still 
sizeable in the considered temperature range; note, how- 
ever, that the simplistic formula (0| does not correctly 
capture these corrections. The small off-set between the 
experimental data and the fit in Fig. |4jd) amounts to a 
systematic error in the effective prefactor 7 at criticality 
of -10% 0. 

Fig. [5] compares the uniaxial magnetostriction coeffi- 
cient A along the c* axis with the ones along the or- 
thogonal a and b directions. While the former is nearly 
symmetric with respect to (H c \ + H c i)j1 this symmetry 
is missing in the other two directions. This difference is 
due to the additional contribution arising from the pres- 
sure dependence of the leg coupling J|| . It turns out that 
the dependence of J\\ on uniaxial pressure p c * \\c* is very 
weak, |9J||/i9p c *| <C \dJ±/dp c *\, such that its influence 
can be neglected. The single dominant p c * -dependent en- 
ergy scale J± causes peaks in A at H cl / c2 that are of com- 
parable sizes as dH cl /dp c * ~ dH c2 /dp c * ex dJ±/dp c * 
[12J ] . This is not the case for the orthogonal directions 
where the respective uniaxial pressure dependencies of 
J\_ and J|| are of similar magnitudes. In particular, the 
small relative size of the peaks at H c \ rj (J± — J\\)/gnB 
can be understood as a partial cancellation of the two 
pressure dependencies dJ\\/dp n — dJ±/dp n yielding a 
small dH c \/ dp n , with n = a,b. Taking the pressure de- 



pendencies of both, J± and J\\, into account, we find 
excellent quantitative agreement with Quantum Monte 
Carlo simulations of the Hamiltonian (fT]), solid lines in 
Fig. O whose details will be presented elsewhere 22[ . 

In summary, we measured the magnetostriction and 
thermal expansion of (CsH^N^CuBr^ For both phys- 
ical quantities we find excellent agreement with calcu- 
lations based on a two-leg spin ladder Hamiltonian (|TJ) . 
The thermal expansion a is critically enhanced as 1/vT 
close to the two quantum critical points at H c \lc2 an d 
shows various sign changes as a function of H and T 
signaling entropy extrema in the phase diagram. This 
complex behavior of a is semi-quantitatively explained 
within a model of free fermions, and we find quantitative 
agreement with Quantum Monte Carlo calculations. 
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